High-frequency electron paramagnetic resonance investigation of the Fe 3+ 
center in polycrystalline PbTiO,3 in its ferroelectric phase 

Hrvoje Mestric, Riidiger- A. EichelE and K.-P. Dinsc 
Eduard-Zintl- Institute, Darmstadt University of Technology, D-64287 Darmstadt, Germany 



impurity 



o 
o 

> 
O 

00 



o 

CO 



Andrew Ozarowski, Johan van Tol, and Louis Claude Brunei 
Center for Interdisciplinary Magnetic Resonance, 
National High Magnetic Field Laboratory, Florida State University, Tallahassee, FL 32310 

(Dated: 2nd February 2008) 

The intrinsic iron(III) impurity center in polycrystalline lead titanate (PbTiOa) was investigated 
by means of high-frequency electron paramagnetic resonance (EPR) spectroscopy in order to deter- 
mine the local-environment sensitive fine structure (FS) parameter D. At a spectrometer frequency 
of 190 GHz, spectral analysis of a powder sample was unambiguously possible. The observed mean 
value D = +35.28 GHz can be rationalized if Fe 3+ ions substitute for Ti 4+ at the B-site of the 
perovskite ABO3 lattice forming a directly coordinated Fe^i — V" defect associate. A consistent fit 
of the multi-frequency data necessitated use of a distribution of D values with a variance of about 1 
GHz. This statistical distribution of values is probably related to more distant defects and vacancies. 



PACS numbers: 61.72.Ji, 61.72.Hh, 61.72.Ww, 76.30.Fc, 77.84.Dy 
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I. INTRODUCTION 

Lead titanate (PbTi03, PT) is widely used as func- 
tional ceramic because of its excellent physical and elec- 
tromechanical propertiesii2i2i4. It is a 'displacive-type' 
ferroelectric material that can be used as dielectric in 
capacitors and as high refractive index thin film for 
electrooptical components. Its piezoelectric and pyro- 
electric properties can be utilized in sensors, piezoelec- 
tric actuators and detectors for infrared radiation. As 
compared to the solid solution system lead zirconate 
titanate(Pb[Zr a; Tii_ x ]03, PZT), lead titanate exhibits a 
higher Curie temperature (Tq = 763 K) and a lower di- 
electric constant of about 200, rendering it more attrac- 
tive for high-temperature and high-frequency transducer 
applications. 

In order to improve material properties, several tran- 
sition metals or rare-earth elements may be added on a 
percentage level, for which reason considerable interest 
exists to characterize the role of such extrinsic functional 
centers. Generally, standard bulk characterization tech- 
niques fail due to the inherent low concentration of these 
centers and electron paramagnetic resonance (EPR) be- 
comes the method- of- choice because of its high sensitiv- 
ity and selectivity^. In particular, if doping with iron 
is considered, the local symmetry can be explored by 
monitoring the resulting fine structure (FS) interaction, 
which will be modified by the presence of oxygen va- 
cancies (Vq*)£. If charge compensation occurs in the 
nearest-neighbor O 2- ion shell, a large distortion of the 
octahedral symmetry results and because of the short 
distance a large change of the intrinsic Fe 3+ FS tensor 



will be induced^. 

Unless crystalline samples are available, these tensor 
elements in general are not directly accessible at X-band 
(9.4 GHz) frequencies, since the resulting zero-field split- 
ting (ZFS) is much larger than the microwave (mw) quan- 
tum energy. FS values therefore have to be deduced from 
the analysis of second-order line shifts in a fully resolved 
single crystal spectrum. If only powder samples are avail- 
able, high-frequency EPR has to be invoked in order to 
approach high- field conditions, under which an accurate 
determination of the principal value D of the FS interac- 
tion (including sign) via first-order effects is possible. 

Beyond its technical importance, we can use lead ti- 
tanate as a model system related to the considerably 
more complex acceptor-doped PZT system^, because its 
structure being comparatively simple and well defined. 
Furthermore, a prerequisite for an analysis of the EPR 
spectra of solid-solution PZT systems with varying Pb/Zr 
composition is a complete characterization of the spec- 
tra of pure compounds. As stated above, in general it is 
quite difficult to extract large ZFS tensor elements from 
X-band EPR powder spectra. In order to establish high- 
frequency EPR as reliable tool for the investigation of 
technologically relevant polycrystalline compounds, we 
have chosen the ubiquitous Fe 3+ impurity center in lead 
titanate powder as first example, because FS parame- 
ters have been obtained by previous single crystal studies 
which can be used for comparison. 

The first EPR spectrum, attributed to Fe 3+ in PbTi0 3 
led to the discovery of an unusual high FS splitting 9 . 
This observation was related to the strong tetragonal fer- 
roelectric distortion of the crystalline field. At Q-band 
(35 GHz), two different iron centers were observed, one 
of them being assigned to a partially charge-compensated 



- V" defect associate. It was also observed that the 
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FS interaction increases by 20 % when cooling from room 
temperature to 77 K, consistent with a change in lattice 
deformation, the c/a ratio changing by f 9 % in this tem- 
perature range. The spectra allowed also for the deter- 
mination of additional terms attributed to higher-rank 
tensor elements. The temperature dependence of line 
intensities at liquid helium temperatures indicated that 
the sign of D is positive!*. At 70 GHz, additional tran- 
sitions could be observed and allowed to refine the ZFS 
parameters^. In contrast to these findings, in a recent X- 
band EPR study only a single Fe 3+ center was observed, 
for which charge compensation was assumed to occur at 
distant spheres!*. This discrepancy was attributed to 
different synthesis techniques involved in manufacturing 
the crystals. At low temperatures hyperfine structure, 
assigned to interaction with nearby 207 Pb nuclei was ob- 
served. 

The reported ZFS parameters are spread over a broad 
range. Since many attempts for gaining structural in- 
formation about the dopant site, such as the assignment 
to either Fe^ — V" defect associates or not-coordinated 
'free' Fe^; centers with or without off-center shifts of the 
iron ion, are based upon size and sign of the ZFS pa- 
rameters, their precise determination is of considerable 
importance. By taking the reported interval of ZFS pa- 
rameters as input parameters for modelling the struc- 
ture of the iron center, a contradictory variety of pro- 
posed structures results, ranging from 'free' Fe^i centers 
to Fe^j — V" defect associates with off-center shifts to- 
wards or away from the oxygen vacancy. In this work, 
we hence aim for an unambiguous determination of the 
ZFS parameter by application of high-frequency EPR up 
to 190 GHz, thus creating the basis for an accurate mod- 
elling of the structure. In the magnetic field range of 2 to 
8 T used, the electron Zeeman energy will be the domi- 
nant term in the spin Hamiltonian. Resulting quantiza- 
tion approximately along the field direction significantly 
simplifies the EPR spectrum. Accordingly it is possi- 
ble to accurately determine ZFS parameters of Fe 3+ in 
polycrystalline PbTi03 from observed van Hove singular- 
ities, which arc related to canonical orientations of the 
compound. High frequency EPR thus is an indispens- 
able tool for the investigation of technologically relevant 
PZT, which are commonly available only as disordered 
compounds. 



II. EXPERIMENTAL 

9.5 GHz continuous wave (c.w.) EPR measurements 
were performed using a ESP 300E spectrometer (Bruker), 
equipped with a rectangular TE112 resonator. The mag- 
netic field was read out with a NMR gaussmeter (ER 
035M, Bruker) and as a standard field marker polycrys- 
talline DPPH with g = 2.0036 was used for the exact 
determination of the resonance magnetic field values. 
High-frequency EPR measurements were performed at 
the National High Magnetic Field Laboratory (NHMFL) 



Tallahassee! 3 -. The set-up used operates in transmission 
mode and employs oversized cylindrical wave guides. No 
resonator was used. Microwave detection was performed 
with a low-noise, fast response InSb hot-electron bolome- 
ter (QMC Ltd.), operated at liquid- helium temperature. 
Field modulation in the range of 10 to 50 kHz was used 
to obtain 'first-derivative' type EPR spectra. 



III. THEORETICAL DESCRIPTION 

In principle, a discussion of the Fe 3+ center in its 
high-spin form (S = |) should be based upon the most 
general form of the spin Hamiltonian, including fourth- 
rank tensor components!*. However, literature values 
of spin-Hamiltonian parameters for Fe 3+ ions in lead ti- 
tanate single crystals, as well as our high-field EPR re- 
sults justify considerable simplifications. First, since no 
significant line broadening was observed in high magnetic 
fields, the anisotropy of the g-matrix apparently imposes 
a much smaller angular-dependent EPR resonance shifts 
than the inherent peak-to-peak linewidth AB pp of the 
resonance lines and might therefore be neglected. This 
is in agreement with expectation because of a vanishing 
orbital momentum of Fe 3+ ions. As a consequence, the 
electron Zeeman interaction will be represented by an 
isotropic g-value. Second, due to the relatively low nat- 
ural abundance of the only magnetic active isotope 57 Fe, 
iron hyperfine interaction will be omitted for spectrum 
reconstruction. Third, the nuclear Zeeman interaction is 
ignored as well because it does not contribute to EPR 
spectral features in first order. Finally, no influence of 
fourth-rank tensor elements could be experimentally re- 
solved. Hence, the corresponding terms were neglected. 

The resulting simplified effective spin Hamiltonian is 
given by 

W = S.D.S + «? iso /3 e B -S (1) 

in which gi SOl g n are the electron and nuclear g factors and 
(3 ei [3 n are the Bohr and nuclear magnetons. The second 
and third terms are the electronic and nuclear Zeeman 
interactions, respectively, and Bo is denoting the exter- 
nal field. The first term describes the second-rank fine 
structure interaction. This term is often referred to as the 
zero-field splitting term because it lifts the degeneracy of 
energy states independent of magnetic field. The FS ten- 
sor D is symmetric and traceless. Thus, there is always 
a coordinate system in which the tensor is diagonal with 
the elements Dm and D±. By convention, D\\ is taken to 
be the principal value with the largest absolute magni- 
tude. Hence, for axial symmetry, the fine-structure split- 
ting may be described in terms of a single crystal-field 
parameter that is commonly defined as D = |Z?|| = 3D±. 
The free Fe 3+ ion has a 3d <S*5/2 ground state. Under 
the action of a tetragonal crystalline field, this ground 
state splits into three twofold degenerate levels (cf. fig- 
ure H(a)). 
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In case of the commonly employed X-band frequen- 
cies, the zero-field splitting is much larger than the elec- 
tron Zeeman energy. All doublets are well separated and 
only transitions within each doublet can be observed. 
Accordingly, the spectrum then is described by defin- 
ing a separate effective spin Hamiltonian with S' = h 
for each doublet. Because of orientation-dependent state 
mixing, the effective ^'-matrices can exhibit a very large 
anisotropy even retaining the assumption that the pure 
electron Zeeman interaction is isotropic 1 ^. This situa- 
tion is schematically represented in figure ^a, c). At 
zero field, the spins are quantized along the z-axis form- 
ing three groups of doubly degenerate states. After ap- 
plication of a static field along the z-axis, the quantiza- 
tion along the crystal z-axis is retained, the degeneracy 
being lifted proportional to Bq. However, due to the 
zero-field splitting being larger than the mw quantum of 
energy at X-band frequencies, only transitions between 
levels degenerate at zero field may be induced, providing 
such transitions are not 'forbidden' by the selection rule 
Amg = ±1. If the magnetic field is perpendicular to the 
z-axis, there is no first-order splitting of the \ms = ±§) 
and | ± |) states, whereas the | ± \ ) states are split with 
g' « 6. Due to the selection rule all transitions, except 
the |ms = +i}^| — transition are 'forbidden' and 
will not be observed in the X-band EPR spectrum. 

IV. RESULTS AND DISCUSSION 

A typical X-band EPR spectrum is presented in fig- 
ure ^b). It mainly consists of one prominent feature at 
low field. In terms of effective g-values, the spectrum 
is described with gj| = 2.004 and g' ± = 5.956. Because 
all orientations are statistically realized in a powder, the 
resonances extend from g',, to g'j_ , as illustrated in figure^ 
(c) , and because the probability of Bo being perpendic- 
ular to the z-axis is largest, the resonance peak at g' ± is 
the dominant feature. The additional feature at g = 2.2 
can be assigned to Pb 3+ centers^. 

The W-band EPR spectrum (95 GHz), depicted in 
figure EJa), is representative for an intermediate field 
regime. For this situation the crystal-field terms are com- 
parable to the electron Zeeman interaction. Hence, no 
simple prediction pattern of the EPR spectrum results. 
The observed sharp lines, however, can be identified as 
arising from to resonances with minor orientation depen- 
dence, as shown in figure Gib). To obtain good agree- 
ment between simulation^ and experiment, it had to be 
assumed that the sample was exposed to oscillating Bi 
field components being both parallel and perpendicular 
to Bo- This phenomenon is caused by multiple reflec- 
tions and general imperfections (mode impurity) in the 
mw propagation system, in contrast to cavity-based sys- 
tems with a well-defined Bi polarization^. For spectrum 
simulations this effect necessitates consideration of paral- 
lel and perpendicular EPR modes. In the perpendicular 
mode, the excitation and detection mw fields are along 




Figure 1: (a) Schematic representation of energy levels as 
function of magnetic field for the canonical parallel (left) and 
perpendicular (right) orientations relative to the ZFS-tensor 
principal axes for the Fe 3+ impurity center in PbTiC>3. The 
FS interaction is taken to be positive, D > 0. EPR transitions 
for X-band frequencies are marked by vertical lines involving 
'allowed' (solid) and 'forbidden' transitions (dashed), (b) X- 
band EPR spectrum of the Fe 3+ impurity center in PbTiC>3 
at v mw = 9.320 GHz and at 10 K. The line at g = 2.2 marked 
by an asterisk is tentatively assigned to Pb 3+ centers, (c) 
Orientation dependence of the resonance lines for 'allowed' 
and 'forbidden' transitions. 



the laboratory cc-axis, whereas in the parallel mode they 
are along the z-axis, parallel to the external static field. 
When the magnetic field is parallel to a principal axis, 
the EPR Ams = ±2 transition probability drops to zero 
with the modulation field Bi perpendicular to Bo, while 
it is largest with Bi parallel to Bo, enabling the obser- 
vation of 'forbidden' transitions. 

At G-band (190 GHz), EPR is performed almost in 
the high-field regime. The corresponding energy-level di- 
agrams for the canonical orientation perpendicular rela- 
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Figure 2: (a) W-band EPR spectrum of the Fe impurity 
center in PbTi0 3 at u m „ = 95.284 GHz (T = 10 K). Numer- 
ical spectrum simulations involve Bi in perpendicular (top) 
and parallel (bottom) mode, (b) Orientation dependence of 
EPR transitions. 



tive to the ZFS-tensor principal axes is depicted in fig- 
ure a). The field-energy dependencies are quite com- 
plicated because of extensive level mixing, resulting in 
strongly varying resonance fields for most of the transi- 
tions. The EPR spectrum is presented in figurc[!Jb). At 
low temperatures, the sign of the D can in principle be 
determined. This is due to the fact that the transition 
intensities are a function of the Boltzmann populations 
of the levels involved. In our case FS splitting and elec- 
tron Zeeman energies are much larger than fc#T at 10 K 
at which temperature spectra were recorded for the de- 
termination of the sign of D. Under these conditions, the 
temperature dependence of particular resonances allows 
to determine the sign of the splitting parameter D. In 
figure |3| (a) , results for numerical spectrum simulations 
assuming positive (solid) and negative sign (dashed) of 
D are shown. Because the resonance at 4.61 T (marked 
by an asterisk) occurs in the calculated spectrum only 
for D > 0, the ZFS interaction can be taken as posi- 
tive. The disappearance of the signal for D < can be 
understood by taking into account the calculated energy 
level diagrams in figure [3^a), where the particular line 
corresponds to the transition connecting states of lowest 
energy. If the sign of D would be reversed, this level 
would be highest in energy and hence this transition has 
to be thermally activated and should be detectable at 
elevated temperatures only. 

The refinement of the spin-Hamiltonian parameters 
through numerical spectrum simulation simultaneously 
for all mw frequencies led to the results summarized in 




(i) 



(iii) 



■A \- 



(iv) 



(v) 



e„[T] 




e [T] 



Figure 3: (a) G-band EPR spectrum of the Fe 3+ impurity 
center in PbTi0 3 at i/ mw = 189.962 GHz and T = 10 K (iii). 
Numerical spectrum simulations involve both Bi in perpen- 
dicular (i,ii) and parallel (iv,v) mode, as well as a positive 
(solid) and negative sign (dashed) of the axial ZFS parameter. 
The transition being particularly sensitive to temperature- 
dependent level population is marked by an asterisk, (b) En- 
ergy levels as function of magnetic field for the canonical par- 
allel (left) and perpendicular (right) orientations relative to 
the ZFS-tensor principal axes using D > 0. EPR transitions 
for G-band frequencies are marked by vertical lines involving 
'allowed' (solid) and 'forbidden' transitions (dashed). 



table I. Using these best-fit values, numerically simu- 
lated spectra are superimposed to the experimental data 
and show excellent agreement. Concerning the relatively 
large value of D, the local environment for the Fe 3+ ion 
can be modelled with the help of the semi-empirical New- 
man superposition mode&^. The geometry of the next- 
nearest oxygen positions was calculated using the lattice 
parameters c = 0.422 nm and a = 0.386 nm from XAFS 
measurements at 12 K — . The parameters for the iron- 
oxygen ion pair were taken as 62 = —12.36 GHz at a 
reference distance of i?o = 0.2101 nm, and the power-law 
exponent as £2 = 82i, respectively. Only by assuming 
Fe 3+ being located at the B-site correlated with a next- 
nearest oxygen vacancy, the measured value for the ZFS 
could be reproduced within the centered model^. An 
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Table I: Table I: g-values and FS parameters as obtained from high-field EPR, compared to values from literature with ^ mw - 
The sign of the principal ZFS parameter is determined by an analysis of line intensities at low temperatures. (For a definition 
of the fouth-order parameters a, F see for example: A. Abragam, B. Bleaney: Electron Paramagnetic Resonance of Transition 
Ions, Clarendon Press, Oxford (1970)) 



center 




D [GHz] 


a [MHz] 


F [MHz] 


T [K] 


!/ mw [GHz] 




Fe' - Vo" 


2.002 


35.28 ±0.01 






2-10 


9.6 - 190 


this work 


Fe' 


2.009(5) 


27.13 ±0.01 
32 ± 1 


1678 ± 10 


-2833 ± 10 


290 
77 


9.4 
9.4 


V2 


Fe' - Vo" 

Fe' 




35.59 ±0.06 
34.48 ± 0.09 


600 ± 210 
660 ± 300 


-510 ± 300 
-600 ± 300 


85 


70 




Fe' - Vo" 

Fe' 

Fe' - Vo" 

Fe' 


2.0 
2.0 


26.98 ± 1.5 
15.89 ±0.6 
32.98 ± 1.5 
20.08 ± 0.6 


809 ± 120 




290 
77 


37.5 


ill 


Fe' 


2.009(5) 


> 30 








9.4 





off-center displacement of the iron in an intact oxygen oc- 
tahedron leads to a considerably smaller value (D < 25 
GHz). Based upon the experimentally determined FS 
value, the conclusion is thus drawn that iron builds a 
defect associate with a directly coordinated oxygen va- 
cancy (Fe^i — V"). These results are in agreement with 
predictions obtained by ab-initio calculations, which are 
in progress. 

Moreover, it can be concluded from the axial symmetry 
of the ZFS tensor that the orientation of the Fe^ ; — V" 
defect dipole is along the crystallographic c-axis. Any 
other coordination of the oxygen vacancy along either a 
or b would reduce the symmetry of the ZFS tensor to 
orthorhombic, which was not observed within spectral 
resolution. 

Starting with D as a single fit parameter, a satisfactory 
reproduction of the experimental data taken at all differ- 
ent mw frequencies was possible only when allowing for 
statistical Gaussian distributions SD. Considering the 
inhomogeneous charge-compensation mechanisms, pro- 
vided by either lead vacancies (Vp h ) or valency-altered 
Pb 3+ centers, and considering that these defects can be 
in different shells of the Fe^ ; — Vq* defect associate, a sta- 
tistical distribution (strain) of the spin-Hamiltonian pa- 
rameters defined by corresponding variances is expected. 
In our case, a variance 8D — 1 GHz of the ZFS parameter 
D, taken to be Gaussian and independent over the mw 
frequency range from 9.6 to 94 GHz, had to be assumed. 

The high-field approach thus allows for determining 



absolute value and sign of the axial fine structure pa- 
rameter for polycrystalline compounds without ambigu- 
ity and even superior accuracy as compared to single 
crystal studies at X-band frequencies. This is due to the 
fact that FS values at high fields can accurately be deter- 
mined via first-order effects rather than by second-order 
shifts at low fields. A second consequence is that a distri- 
bution 8D of the FS parameter is observable as first-order 
variation of the line positions that result in line broad- 
ening in the high-field spectra. At X-band, this effect 
is of second-order only and thus almost vanishes. Hence 
no information about £)-strain can be gathered. In con- 
trast, fourth order FS parameters were not accessible by 
numerical spectrum simulations up to 190 GHz, which 
in turn are accessible by orientation-dependent single- 
crystal studies at X-ban d 10 ! 11 ' 12 . 
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